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Transforming growth factor (TGF) /3 is known to have properties of both a tumour 
suppressor and a tumour promoter. While it inhibits cell proliferation, it also increases cell 
motility and decreases cell-cell adhesion. Coupling mathematical modeling and experiments, 
we investigate the growth and motility of oncogene-expressing human mammary epithelial 
cells under exposure to TGF— /3. We use a version of the well-known Fisher-Kolmogorov 
equation, and prescribe a procedure for its parametrisation. We quantify the simultaneous 
effects of TGF-/3 to increase the tendency of individual cells and cell clusters to move 
randomly and to decrease overall population growth. We demonstrate that in experiments 
with TGF-/3 treated cells in vitro, TGF-/3 increases cell motility by a factor of 2 and 
decreases cell proliferation by a factor of 1/2 in comparison with untreated cells. 

Keywords: Transforming growth factor (TGF) /3, cell growth and motility, mathematical 
model, Fisher-Kolmogorov equation 
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1. Introduction to the Biology of TGF /3 

In normal organisms, the growth of cells is under tight regulation by growth fac- 
tors and is highly dependent on the developmental stage during the lifespan of the 
organism. Disruption of this regulation is the most frequent cause of cancer dis- 
eases. Unlike normal differentiated cells, cancer cells are usually hyperproliferative 
as a result of the abnormal activation of multiple growth-stimulating intracellular 
signalling pathways and loss of tumour suppressors. In cancer cells, these pathways 
do not respond to normal regulatory signals but are manipulated by one or more 
oncogenic signals, often encoded by oncogenes. Expression of oncogenes alters sig- 
nalling pathways that under normal conditions maintain cell growth homeostasis. 
Thus, altering these pathways may favour increased cell and cancer growth. A good 
example is the transforming growth factor (TGF) (3 family, which is known to be 
able to act as both a tumour suppressor and tumour promoting factor. 

The TGF-/3 family consists of multitasking cytokines that play important roles 
in cell proliferation, cell motility, apoptosis, lineage determination, extracellu- 
lar matrix production, and modulation of immune function [22 . These ligands 
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bind to a heteromeric complex of transmembrane serine/threonine kinases, the 
type I and type II receptors (T/3RI and T/3RII). The receptors are activated 
upon hgand binding, leading to the subsequent phosphorylation and activation 
of a family of transcription factors called Smads, which regulate transcription 
of a subset of genes |23j. In addition to Smads, other signalling pathways have 
been implicated in TGF-/3 actions in recent studies. These include the extracel- 
lular signal-regulated kinase (ERK, MAPK), c-Jun NH2-terminal kinase (JNK), 
p38MAPK, phosphatidylinositol-3 kinase (PI3K), and Rho GTPases (reviewed in 
[71 [TTl |39l Us]). The critical role of these non-Smad pathways on mediating the 
cellular effects of TGF-/3 remains to be fully characterised. 

TGF-/3 was originally reported to induce transformation of mouse fibroblasts 
[23]. Subsequent studies indicated that TGF-/3 is a potent inhibitor of cell pro- 
liferation and a tumour suppressor |321 136] . Consistent with its tumour suppres- 
sor role, many cancers lose or attenuate TGF-/3-mediated anti-mitogenic action 
by mutational inactivation of TGF~/3 receptors or their signal transducer Smads 
[T3l [T5l [T6l l20l liOl W\\ . There is increasing evidence to show that excess production 
and/or activation of TGF-/3 in tumours can accelerate cancer progression through 
enhancement of tumour cell motility and survival, increase in tumour angiogenesis, 
extracellular matrix production and peritumoural proteases, and the inhibition of 
immune surveillance mechanisms in the cancer host (reviewed in [3 |9l |Tl] ) . Cancer 
progression and metastasis consist of a series of sequential events. After initial cell 
transformation, often mediated by the function of oncogenes, tumour cells grow- 
ing at the primary site will invade the surrounding stroma and migrate towards 
blood vessels. Through various mechanisms such as epithelial-mesenchymal tran- 
sition (EMT), tumour cells will enter the blood vessels and travel to other parts 
of the body through the circulatory system. Some of the cells will then arrest at 
distant sites where they may proliferate and invade into the adjacent organs. Cell 
motility is therefore a critical element during the spread of tumour cells from their 
initial sites of residence. In this study, we focus on the tumour-promoting effect of 
TGF-/3 through inducing cell motility. 

The receptor tyrosine kinase HER2 (ErbB2, Neu) belongs to the family of epi- 
dermal growth factor receptor (EGFR). Gene amplification or overexpression of 
IIER2 is observed in about 25% of breast cancers. TGF-/? has been shown to 
synergize with the oncogene ErbB2 in cancer progression. Overexpression of ac- 
tive TGF-/3 1 or active mutants T/3RI (Alk5) in the mammary gland of bigenic 
mice also expressing mouse mammary tumour virus (MMTV)/Neu (ErbB2) accel- 
erates metastases from Neu-induced mammary cancers |25] [26l |271 135] . Exogenous 
as well as transduced TGF-/3 confer motility and invasiveness to MCFIOA non- 
transformed human mammary epithelial cells (HMEC) stably expressing trans- 
fected IIER2 |34| [38] . Expression of the oncogene HER2 in these cells does not 
affect the function of TGF-/3 on inhibiting cell proliferation ^38j. It is likely that 
in many cancers, TGF-/3 may still attenuate proliferation while inducing cellular 
events associated with metastatic dissemination, such as cell motility. In this paper 
we report experiments with MCF10A/IIER2 cells to study and to separate the ef- 
fects of TGF-/3 on cell proliferation and motility. Due to the complexity of TGF-/3 
signalling that simultaneously affects several biological parameters, it is important 
to computationally simulate the behaviour of cells under TGF-/3 exposure. Our 
model, which can also be adopted to simulate other growth-regulating signals, will 
provide a unique insight into the TGF-/? function in both normal and cancer cells, 
and further understanding of targeted therapeutic strategies that aim at interfering 
with TGF-/3 signalling. 

Mathematical modeling of chemotherapy with Tamoxifen has been carried out 
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in [37], with special consideration of TGF-/3. There, the authors use a discrete 
cellular automaton (CA) model to investigate the effects of TGF-/3 on tumour 
morphology and invasiveness. It is possible to assign properties to individual cells 
in a cellular automaton model (such as different levels of malignancy). However, 
they can be computationally costly, in particular in situations where large cell 
populations are involved. Since we are interested chiefly in the effects of TGF-/3 on 
a large number of cells growing in a plate, rather than individual cell behaviour, 
we work with a continuous partial differential equation model. Future research 
problems may require the construction of hybrid models such as in [Ij where discrete 
cells are coupled to continuous fields that in turn are solutions of partial differential 
equations. 



2. Introduction to the Mathematical Model 

The Fisher-Kolmogorov equation (or KPP equation) has been studied and used 
widely in mathematical biology. It originated in the 1930s in works of R. A. Fisher 
and A. N. Kolmogorov |12l [T8] where it was proposed as a model for the spread of 
an advantageous gene in a population. Mathematically, it falls into the large class 
of reaction-diffusion equations, in which one or more diffusible "species" enter 
into a scheme of reactions. The species can be chemical substances or biological 
populations. Reaction-diffusion equations in general and the Fisher-Kolmogorov 
equation in particular have been applied successfully by many authors to model 
the behaviour of cells in tissues, see for example [U 13 |3 El El EHl ESI ED]- The 
Fisher-Kolmogorov equation is usually used to show traveling wave behaviour, but 
we use it here to model the spatial growth of a cell culture from initial seeding to 
confluence. 

Let ii(x, y, t) denote the spatial density of tumour cells at time t in a dish corre- 
sponding to a spatial region 17 C K^. The time evolution of this density is described 
by the partial differential equation 



Here a incorporates the intrinsic proliferation rate (unit time~^), from initial seed- 
ing to confluence, in a single value independent of x, y, t. Also, /? denotes the area 
occupied by an average single cell. The density satisfies, u{x,y,t) < (3^^ for all 
X, y, t, and the total number of cells in the region 17 C at time t is given by 



The spatial region $7 may be the entire region occupied by the population or a 
subfield thereof. Equation ([T]) is accompanied by an initial condition and appro- 
priate boundary conditions. If a subfield of O is completely populated, then 
u{x,y,t) = for all x,y E Oi, otherwise u{x,y,t) < I3~^ and 



d_ 
di 



u{x, y, t) = DAu{x, y, t) + au{x, y, t){l — f3u{x, y, t)). 



(1) 





is the fraction of possible cells in at time t. Here and in the sequel | • | denotes 
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the area of the domain in question. We denote the total mass of the density by 

U{t)=p [ [ u{x,y,t)dxdy. (2) 
J Jn 

Much of the remainder of this paper wih be concerned with the problem of finding 
parameter values to a modified version of equation ([T| from experimental observa- 
tions and with the interpretation of the solution to the equation. 

The diffusion constant D in equation ([T]) has to account for three simultaneous 
processes contributing to the spatial movement of cells. First, dividing cells occupy 
increasing space, therefore a certain part of D, say Dp, has to account for the 
spatial expansion of proliferating cells. The authors of [6l 13 [191 EO] have proposed 
relations of the type 



Dp = ki 



Td 



Here i is related to the typical length or diameter of a cell and is the doubling 
time. The dimensionless factor ki in [S^ accounts for the expansion of the viable 
rim of dividing cells. We take (. to be the average increase in cell diameter between 
mitotic events. We assume the area of an average cell increases from approximately 
50 ixm^ to WOfim?, so we take i = 3.3 nm. We assume the average cycle time is 
Td = 16 h. We take Dp = 0.1 iJ,m?h~^ (order of magnitude) . Second, individual cells 
plated on a Petri dish do not remain fixed in their position but undergo a random 
walk. It is natural to assume that the random walk is not biased and follows the 
laws of Brownian motion. Third, cells in a cluster may break lose from that cluster 
as TGF-/3 is known to decrease cell-cell adhesion |28j. Thus, a part of D, say Dm, 
must account for this second and third contribution to cell movement. 

We propose here a modification of the standard Fisher-Kolmogorov equation 
^. Our reasoning is that as the cells become more densely packed, their random 
motility should decrease. Thus, with U from equation ([2| we propose 

d 

Q^u{x, y,t) = V ■ {D{U)Vu{x, y, t)) + au{x, y, t){l - /3u{x, y, t)), 
D{U) = Dm[l-^^j+ Dp. 

The constant Dm, which captures the random motility of individual cells and clus- 
ters and the tendency of clusters to break apart, depends on TGF-/3 concentration. 
We remark here that D could also be made dependent on the local density u. In 
the case of interest to us in the present paper the initial seeding of the cells results 
in a uniform dispersion of the population that justifies the approximation made 
in equation ([s]). Since Dp > 0, D{U) > and the modified equation (|3| has the 
important property of mathematical well-posedness. In practical applications. Dp 
is often considerably smaller than Dm (see the Results section and table [2] for our 
numerical parameter values). 

Different numerical values for the constant Dm have been proposed in the liter- 
ature. Bray reports in his 1992 book [4, Table 1-1] Dm = 5 • W~^^ cm^ = 
180 fim'^h^^ . Chaplain and Matzavinos [6] arrive at a similar value Dm = 7 • 
10~ cm" day-^ = 300 ^m?h ^. This value is obtained in [6! from the Einstein- 
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Stokes equation |3] 

Svrar/ 

where ks is Boltzmann's constant, T the temperature, d the diameter of the cells, 
and r] the dynamic viscosity of the surrounding medium. It can be excluded, how- 
ever, that cells in vitro, such as in our experiments, move in an environment that 
has the viscosity of water, tjh^o = Pas. Thermally driven motion and Stokes 
friction are not relevant in the context of moving cells. 

3. Materials and Methods 

MCF10A/HER2 cells were generated and grown as described |38j. All live cell 
imaging experiments were performed in triplicate. For growth assay, equal number 
of cells (1.5 • lO'^/well) suspended in full medium were seeded on 6-well plates. 
Cells were allowed to grow in the absence or presence of TGF-/? at different con- 
centrations (0.5 — bngml^^) over a time course of 8 days. Medium containing 
fresh TGF-/? was replenished every 2 days. Cells were harvested by trypsinization 
every 24 h, and subjected to total cell number counting using a Coulter counter. 
In Figure |4]B each data point represents the mean of 3 wells. Standard deviations 
were less than 5% in all cases and were not included in the figure. For random 
motility assay of individual cells, cells were seeded on 6-well plates at 3 • 10^/well 
one day prior to image recording. For motility of cell clusters, cells were seeded 
at 10^/well two days prior to image recording. Culture media was replaced with 
fresh serum-free or full L-15 medium and different concentrations of TGF-/3 was 
added 1 h before the image recording. Cells were imaged on a Nikon TE 2000-E 
microscope and captured using a Hamamatsu Orca ER camera and Metamorph 
software. Images were taken every 5 minutes for 270 minutes for random motility 
assays and every 30 minutes for 14 hours for the cell cluster experiments. Image 
analysis was performed using Metamorph software (MDS, Inc., Toronto, Canada). 

4. Results 

4.1. Random motility of individual cells 

We have taken position measurements of individual cells and have obtained time 
series data {xi,yi)^i, with AT = brain. The mean-squared displacement (MSD) 
is calculated according to 

^ N-k 

r\k^T) = Yl - x^f + ivi+k - Vi?) , (4) 

i=l 

where k = 1,2, ... ,kmax- The maximal multiple kmax for which is calculated 
has to be chosen carefully. Obviously, the bigger k, the fewer such displacements 
are contained in the time series of length A'^ thus the larger the uncertainty in the 
estimate r^(A;). The paper [31] gives estimates of the standard deviation of r'^{k) 
in terms of N and kmax- To obtain the random motility of a single cell we fit the 
estimate (|4]) with a linear function (see |3]| for a discussion of this relation) 



r'^ikAT) = 4DmkAT. 



(5) 
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Unfortunately, the analysis is made difficult by cells that seem to remain immobile 
or whose r^(A:)-curves do not look linear. For the plots in Figure [l]we have selected 
only curves that could be fit reasonably by a straight line. Other curves show cells 
moving initially and then coming to a standstill. Moreover, we cannot see a clear 
increase of as TGF-/? is added. For the average random motility of individual 
isolated mobile cells we obtain, both for cells treated with and with bngml~^ 
TGF-/3, 

Dm = 120/xm2/i~^ 

This value can be taken as an upper bound which is realised only if all cells are 
mobile. 

We now look at the percentage of mobile cells. We define a cell to be mobile if it 
moves outside of a 100 ^.m x 100 jjim square centred at the cell's original position. 
It becomes clear (Table [T]) that the percentage of mobile cells in our experiment 
increases as TGF-/3 is introduced. 

It seems to be a frequently observed phenomenon that a certain subpopulation of 
cells remains immobile. Selmeczi et al. [3H] observed different types of keratinocytes 
and fibroblasts and selected trajectories only for cells that were moving at all. This 
selection procedure is clearly admissible for single particle tracking experiments, 
but is not applicable to a population with a high fraction of non-mobile cells. We 
want to emphasize that the cells are not subject to a directional gradient of TGF- 
/3, rather a uniform concentration applied from all directions, therefore they are 
able to move randomly in any direction, as opposed to being stimulated to move in 
one particular direction. Although it is unknown why cells suddenly stop moving, 
it is possible that in the absence of a strong directional signal there is "noise" in 
the motility signaling that leads to immotility in randomly moving cells. 



4.2. Random motility of cell clusters 

We next investigated the motility of cell clusters (~ 10 cells) in the absence or 
presence of TGF-/3. As shown in Figure [2] top row, untreated cell clusters are 
relatively immobile and seem to maintain a certain shape, although the cells can 
still rotate within the clusters. In contrast, clusters treated with TGF-/3 are sig- 
nificantly more active in altering their shape and migrating towards unoccupied 
space (Figure |2] bottom row). We analysed the movie data in the following way. 
On the initial frame of the movie the centre of the cluster is identified as the centre 
of the smallest rectangle containing the cluster. Then, on each subsequent frame of 
the movie, a straight line is drawn to the point on the boundary of the cluster the 
farthest from the initial centre. The time gap between two frames is AT = Ih. We 
obtain a star of radii as shown in Figure A. For a time series of N radii {ri)^^-^ 
we calculate the variation of the squares 

-'=J^i:\rU-rl\. (6) 
1=1 

For a total of 36 movies, 18 control cases and 18 treatment cases (5 or W ng ml~^) 
respectively, a difference can be seen in the value of the variation v'^ (Figure [sjB) 
with TGF-/3 treated cells showing a higher value than untreated cells. Movies were 
obtained both under serum-free conditions and with serum. The clusters treated 
with TGF-/3 are more active by a factor of two. 
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We combined the observations of individual mobile cells and clusters of cells to 
obtain values of Dm = 5 fim'^h^^ (TGF-/3 concentration 0) and Dm = 10 ^m?h^^ 
(TGF-/? concentration bngml^^). These values are scaled downwards from the 
upper bound obtained from the single cell motility experiments in view of the fact 
that only a fraction of cells are fully mobile. The ratio of the Dm's for treated and 
untreated case is suggested by our experimental observations. 

4.3. Growth assays and their numerical simulation 

The mathematical model from equation ([s]) was programmed using matlab (ver- 
sion 7.1, The Math Works, Inc., Natick, MA). A standard Crank-Nicolson scheme 
is used for the diffusion step [291 Section 17.3], combined with an Euler forward 
reaction step. The codes will be available upon request from the corresponding 
author. 

We fix the following length and time scales for our simulations. The square 
domain has side length L = 300 and the time unit is T = Ih. We assume 
that a single cell occupies an area of /3 = WOfim?, thus the field has a carrying 
capacity of 900 cells. We impose periodic boundary conditions to account for the 
fact that our numerical domain is a small section of a larger domain, say a Petri 
dish. Thus, a cell mass that leaves the region will be balanced by a cell coming in 
at the opposite side. It should also be pointed out that the carrying capacity is the 
same in control and treatment cases. 

Numerical simulations of cells growing in a field are shown in Figure [4] A. The 
initial datum resembles the random placement of cells on the Petri dish. Using 
our modified Fisher-Kolmogorov equation Q we can simulate the experimentally 
determined growth data (Figure |4]B, discrete symbol represent experimental data, 
continuous curves represent simulations). At first, the experimental growth data 
were normalised with respect to a maximum capacity of 1.1 • 10^ cells. The parame- 
ter a was obtained from fitting the experimental data (see Fig.|4]B), while Dm and 
Dp were chosen as stated above. The numerical simulation starts only 72 h into the 
experiment. In our previous experiments with these cells p3], we have noticed that 
after the cells suspended in medium are seeded onto the plates, they require about 
48 hours to adhere to and spread on the plates and adjust to their new environ- 
ment, before they start a typical growth regime. The numerical parameter values 
are collected in table |2] To support our choice of a we note that if the cells were 
able to grow exponentially according to the law u' = au, then the corresponding 
doubling times T2 = ln2/a would be T2 = lOh for untreated cells (a = 0.07 /i~^) 
and = 17 h for cells treated with TGF-/? (a = OMh'^). These doubling times 
are feasible for MCF10A/IIER2 cells in the early exponential growth phase. 



5. Discussion and Conclusions 

Our experiments demonstrate that TGF-/3 increases the percentage of mobile cells 
in a cell population in a dose-dependent manner, rather than increases the mean 
square migration displacement of individual cells. As we demonstrated, random mi- 
gration of cell clusters containing ~ 10 cells is a feasible approach to parametrise 
the unbiased random cell migration in a large population of cells. The cluster 
motility assay indicates that clusters are more mobile and/or less cohesive if TGF- 
P is present. In fact, biological evidence for both mechanisms has been identified. 
TGF-/3 induces a PI3K~mediated activation of Racl/Pakl pathway, which leads 
to increased cytoskeleton reorganisation, turnover of focal adhesion, and eventu- 
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ally cell migration |32]. Furthermore, TGF-/3 downregulates cell-cell adhesion by 
decreasing the adhesion protein E-cadherin [28^ . The cluster motility assay reflects 
an overall efl'ect of the multifunctional molecule TGF-/3 on motility. In metastatic 
cancer cells, TGF~/3 stimulates epithelial-mesenchymal transition (EMT) which 
is necessary for the acquisition of invasive and metastatic phenotype. It has been 
shown that in cell monolayer wounded with a pipette tip, the presence of TGF-/3 
in the medium markedly induces the closure of wounded area, which is another 
indicator of TGF-/3-induced cell motility. However, this assay is not used in our 
parametrisation as the cell migration is not random but targeted to the wounded 
area. 

Our model has a relatively small number of parameters: Dp, /?, Dm.) and a. 
Our procedure for determining these parameters is as follows: (1) Dp and (3 are 
determined by the average doubling time and average area of cells in the culture, re- 
spectively. (2) Dm is determined by adjusting the motility measurements of isolated 
mobile cells and clusters of mobile cells for their fraction of the total population as 
the culture attains confluence. From our experiments the upper estimate on D^ 
was 120 jim?h~'^, which we adjusted downward by ~ 1 order of magnitude to com- 
pensate for the fraction of mobile cells in the total population. Table [T] provides the 
fractions of mobile cells for varied concentrations of TGF-/3. (3) a is determined 
by fltting model simulations, with the parameters Dp, (5, and Dm as above, to the 



It is possible that different values for Dm and a give rise to very similar total pop- 
ulation growth curves. We determine ratios of Dm and a for treatment and control 
cases that simulate our experimental growth data (Figure |4]B). Our determination 
of Dm (corresponding to motility) and of a (corresponding to proliferation) reveals 
that TGF-/? roughly doubles the value of Dm and halves the value of a. The total 
population growth data of the cell cultures alone (without any further information) 
are insufficient to quantify these two separate effects of TGF-/3. We have used in 
addition the motility experiment data to distinguish the two effects. 

We have developed a general spatial model of proliferating cell cultures in vitro, 
which allows quantification of the properties of proliferative capacity, cell mobility, 
and clustering as the population attains confiuence. The novelty of our interpre- 
tation of the Fisher-Kolmogorov equation is that cells begin as isolated geometric 
regions corresponding to seeding. As they proliferate and migrate these regions ex- 
pand as controlled by the diffusion parameters and proliferation parameter. When 
the density is identically (whose value is readily correlated to the average area 
occupied by an individual cell), the region occupied is completely populated. When 
the density is < /3~^ in a region, then the integral of the density over this region is 
the expected value of the number of cells in the region. Our focus here is on quan- 
tification of the effects of TGF-/3 on motility and proliferation for in vitro cancer 
cell cultures, but our model allows investigation of similar phenomena for many cell 
types, both prokaryotic and eukaryotic. Extensions and modifications of our model 
could incorporate different cell proliferation dynamics, multiple cell types and cell 
interactions (for example endothelial cells, hematopoietic cells, and immune cells), 
in both in vitro and in vivo settings. In the future we plan to add more components 
that naturally exist in the tumour microenvironment to the model, including stro- 
mal fibroblasts, immune cells and vascular endothelial cells, which are all affected 
by TGF-/3. In addition, many tumours carry key mutations on oncogenes or tu- 
mour suppressor genes, which may affect their responses to TGF-/5. For example, 
cancer cells can become unresponsive to the anti-proliferative function of TGF- 
/9, while remaining sensitive to TGF-/3-induced motility. Inactivating mutations 
of TGF-/? receptors are frequently found in pancreatic, biliary and colon cancers 
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|14| [2T| Ho] . Usually, only a subset of cancer cells in a tumour harbour these muta- 
tions and these cells live together and communicate with other wild-type epithelial 
cells, fibroblasts and endothelial cells. It will be important to further investigate 
with the proper assistance from a mathematical model how these cancers are af- 
fected by TGF-/3. The dual function of TGF-/3 as both a tumour suppressor and 
promoting factor is reflected by its individual effects on inhibiting cell proliferation 
and stimulating cell motility. Primary tumour growth by active cell proliferation 
results in the increase of tumour size. However, during cancer progression, dis- 
semination of tumour cells from primary site into circulation and the subsequent 
invasion of other organs, in which cell motility plays a critical role, usually deter- 
mine the grade of malignancy and the outcome of cancer treatment. Therefore, 
despite the fact that TGF-/3 has been considered a promising therapeutic target 
to treat certain types of cancer, the uncertainty of the overall effects of TGF-/3 
intervention due to the potential risk of derepressed tumour cell proliferation has 
been an unsolved critical issue. This requires further in-depth understanding of the 
different aspects of TGF-/3 action through highly quantitative approaches such as 
mathematical modeling. Our model described here is the first step of encapsulat- 
ing the multiple functions of a signalling molecule and potential therapeutic target 
such as TGF-/3 into a computational model. By introducing the modeling method 
to the cancer biology of TGF-/3, questions that have been raised for years will have 
the chance to be answered from a novel angle. These interesting topics include: (1) 
When is the critical time of TGF-/3 being switched from a tumour suppressor to 
a promoter? (2) When is the most efficient time to apply anti-TGF-/3 therapy 
during cancer progression? (3) Will therapeutic inhibition of TGF-/3 at pre-cancer 
lesions or at early stages of the disease trigger primary tumour growth as a result 
of derepression of cell proliferation? Besides the complexity of TGF-/3 signalling, 
cancer cells can respond differently to TGF-/3. For example, several oncogenes in- 
cluding HER2 are known to synergize with TGF-/3 to function on cell motility and 
survival. Some cancers do not respond to the anti-proliferative effect of TGF-/3 
but still respond to its stimulatory effect on motility. The mathematical model 
will therefore provide a unique way to simulate how TGF-/3 functions differently 
in cancer cell populations with different characteristics, and the effect of TGF-/3 
therapeutic inhibition on a tumour consisting of heterogeneous cell populations. 

6. Acknowledgements 

This work was supported by the NCI Integrative Cancer Biology Program (U54 
CA113007), K22 CA109590 (AMW), K99 CA125892 (SEW), ROl CA80195 (CLA), 
ROl CA62212 (CLA), Breast Cancer Specialised Program of Research Excellence 
(SPORE) P50 CA98131, and the Vanderbilt Integrative Cancer Biology Center 
(VICBC). The authors thank Dr. Vito Quaranta for support. 



References 

[1] Anderson, A. R. A., A hybrid mathematical model of solid tumour invasion: the importance of cell 

adhesion, IMA Math. App. Med. Biol. 22 (2005), pp. 163-186. 
[2] Ayati, B. P., G.F. Webb, and A.R.A. Anderson, Computational methods and results for structured 

multiscale models of tumor invasion, Multiscale Model. Simul. 5 (2006), pp. 1-20. 
[3] Berg, H. C, Random Walks in Biology, Princeton University Press, Princeton, NJ (1983). 
[4] Bray, D., Cell Movements, Garland Publishing, New York (1992). 

[5] Chaplain, M. A. J. and A.R.A. Anderson, Continuous and discrete mathematical models of tumor- 
induced angiogenesis. Bull. Math. Biol. 60 (1998), pp. 857-899. 

[6] Chaplain, M. A. J. and A. Matzavinos, Chaplain, M. A. J. and A. Matzavinos, A. Friedman, Lecture 
Notes in Mathematics, vol. 1872, Springer Verlag, Berlin, New York. 



40 Computational and Mathematical Methods in Medicine TGFbeta'May'l 



10 Taylor & Francis and I. T. Consultant 



[7] Derynck R., R.J. Akhurst, and A.J. Balmain, TGF~f) signaling in tumor suppression and cancer 

progression, Nat. Genet. 29 (2001), pp. 117-129. 
[8] Drasdo, D. and S. Hohme, Individual-based approaches to birth and death in avascular tumors, Math. 

Comput. Model. 37 (2003), pp. 1163-1175. 
[9] Dumont, N. and C.L. Arteaga, Targeting the TGF (3 signaling network in human neoplasia, Cancer 

Cell 3 (2003), pp. 531-536. 
[10] Endcrling, H. ct al., A mathematical model of breast cancer development, local treatment and recur- 
rence, J. Thcor. Biol. 246 (2007), pp. 245-259. 
[11] Ewan, K. B. ct al., Latent transforming growth factor-fi activation in mammary gland: regulation by 

ovarian hormones affects ductal and alveolar proliferation. Am. J. Pathol. 160 (2002), pp. 2081-2093. 
[12] Fisher, R. A., The wave of advance of advantageous genes, Ann. Eugenics 7 (1937), pp. 355—369. 
[13] Gobbi, H. et al., Genetic alterations of the transforming growth factor receptor genes in pancreatic 

and biliary adenocarcinomas. Cancer Res. 58 (1998), pp. 5329-5332. 
[14] Goggins, M. et al., Genetic altemtions of the transforming growth factor /3 receptor genes in pancreatic 

and biliary adenocarcinomas. Cancer Res. 58 (1998), pp. 5329—5332. 
[15] , Transforming growth factor (3 and breast cancer risk in women with mammary epithelial 

hyperplasia, J. Natl. Cancer Inst. 91 (1999), pp. 2096-2101. 
[16] Hahn, S. A. et al., Homozygous deletion map at 18q21.1 in pancreatic cancer, Cancer Res. 56 (1996), 

pp. 490-494. 

[17] Hinow, P. et al., A rnaihernatical model separates quantitatively the cytostatic and cytotoxic effects 

of a HER2 tyrosine kinase inhibitor, Thcor. Biol. Med. Model. 4 (2007), p. 14. 
[18] Kolmogorov, A. N., I.G. Pctrovskii, and N.S. Piskunov, A study of the equation of diffusion with 

increase in the quantity of matter, and its application to a biological problem. Bull. Moscow Univ. 

1 (1937), pp. 1-25, english translation in P. Pelce, editor, Dynamics of Curved Fronts, San Diego, 

1988. Academic Press. 

[19] Lefever, R. and T. Erneaux, Lefever, R. and T. Erneaux, W.R. Ross and A.F. Lawrence, Plenum 

Press, New York. 

[20] Markowitz, S. ct al., Inactivation of the type II TGF-fi receptor in colon cancer cells with microsatel- 

lite instability. Science 268 (1995), pp. 1336-1338. 
[21] , Inactivation of the type II TGF-f3 receptor in colon cancer cells with microsatellite instability, 

Science 268 (1995), pp. 1336-1338. 
[22] Massaguc, J., TGF-f3 signal transduction, Annu. Rev. Biochem. 67 (1998), pp. 753-791. 
[23] Massaguc, J. and Y.G. Chen, Controlling TGF-f} signaling, Genes Dev. 14 (2000), pp. 627-644. 
[24] Moses, H. L. ct al.. Transforming growth factor production by chemically transformed cells, Cancer 

Res. 41 (1981), pp. 2842-2848. 
[25] Muraoka Cook, R,. S., N. Dumont, and C.L. Arteaga, Dual role of transforming growth factor (3 in 

mammary tumorigenesis and metastatic progression, Clin. Cancer Res. 11 (2005), pp. 937s— 943s. 
[26] Muraoka-Cook, R. S. ct al.. Activated type i TGF f3 receptor kinase enhances the survival of mammary 

epithelial cells and accelerates tumor progression. Oncogene 25 (2006), pp. 3408—3423. 
[27] Muraoka, R. S. ct al.. Increased malignancy of Neu-induced mammary tumors overexpressing active 

transforming growth factor 1, Mol. Cell. Biol. 23 (2003), pp. 8691-8703. 
[28] Nawsliad, A. ct al.. Transforming growth factor-beta signaling during epithelial-mesenchymal trans- 
formation: implications for embryoqenesis and tumor metastasis. Cells Tissues Organs. 179 (2005), 

pp. 11-23. 

[29] Press, W. H. ct al.. Numerical Recipes in C, Cambridge University Press, Cambridge (1988). 

[30] Prigogine, I. and R. Lefever, Stability problems in cancer growth and nucleation, Comp. Biochem. 

Physiol. B 67 (1980), pp. 389-393. 
[31] Qian, H., M.P. Shectz, and E.L. Elson, Single particle tracking. Analysis of diffusion and flow in 

two-dimensional systems, Biophys. J. 60 (1991), pp. 910-921. 
[32] Roberts, A. B. et al., Type 3 transforming growth factor: a bifunctional regulator of cellular growth, 

Proc. Natl. Acad. Sci. USA 82 (1985), pp. 119-123. 
[33] Selmcczi, D. ct al.. Cell motility as persistent random motion: theories from experiments, Biophys. 

J. 89 (2005), pp. 912-931. 

[34] Seton-Rogcrs, S. E. et al.. Cooperation of the ErbB2 receptor and transforming growth factor (3 in 
induction of migration and invasion in mammary epithelial cells, Proc. Natl. Acad. Sci. USA 101 
(2004), pp. 1257-1262. 

[35] Sicgel, P. M. and J. Massaguc, Cytostatic and apoptotic actions of TGF-13 in homeostasis and cancer, 
Nat. Rev. Cancer 3 (2003), pp. 807-821. 

[36] Tucker, R. F. et al., Growth inhibitor from BSC-1 cells closely related to platelet type beta transform- 
ing growth factor. Science 226 (1984), pp. 705-707. 

[37] Turner, S., J. A. Sherratt, and D. Cameron, Tamoxifen treatment failure in cancer and the nonlinear 
dynamics of TGF-(3, J. Thcor. Biol. 229 (2004), pp. 101-111. 

[38] Ucda, Y. ct al., Overexpression of IIER2 (erbB2) in human breast epithelial cells unmasks trans- 
forming growth factor (3 -induced cell motility, J. Biol. Chem. 279 (2004), pp. 24505-24513. 

[39] Wakefield, L. M. and A.B. Roberts, TGF-(3 signaling: positive and negative effects on tumorigenesis, 
Curr. Opin. Genet Dev. 12 (2002), pp. 22-29. 

[40] Wang, D. ct al., Analysis of specific gene mutations in the transforming growth factor-(3 signal 
transduction pathway in human ovarian cancer. Cancer Res. 60 (2000), pp. 4507-4512. 

[41] Wang, J. et al., Demonstration that mutation of the type II transforming growth factor (3 receptor 
inactivates its tumor suppressor activity in replication error-positive colon carcinoma cells, J. Biol. 
Chem. 270 (1995), pp. 22044-22049. 

[42] Wang, S. E. ct al., HER2/Neu (ErbB2) signaling to racl-pakl is temporally and spatially modulated 
by transforming growth factor (3, Cancer Res. 66 (2006), pp. 9591-9600. 

[43] Yang, Y. A. et al.. Lifetime exposure to a soluble TGF-13 antagonist protects mice against metastasis 
without adverse side effects, J. Clin. Invest. 109 (2002), pp. 1607-1615. 



May 1, 2008 15:40 Computational and Mathematical Methods in Medieine TGFbeta'May'l 



Computational and Mathematical Methods in Medicine 11 
Table 1. The percentage of mobile cells in the single cell experiments as function of TGF— /3 concentration 
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Tabic 2. The parameter values used in the simulations shown in figure [4| B 
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Figure 1. A Selected vs. k curves for untreated cells. Position measurements of cells were taken and 
the mean-squared displacement was calculated according to equation (j4j. Out of 20 available curves six 
showed a roughly linear dependence of r-^(fe). B Selected vs. k curves for cells treated with TGF-/3. 
The time unit is 5 min for both panels, while the mean-squared displacement is in fim^. 




Figure 2. Random motility of cell clusters in the absence and presence of TGF-/3. The movies are 
represented by 5 different frames at the indicated times. The bar equals 10 fim. 
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A X (n m) B # of movie 

Figure 3. A Coverage radii for cell clusters in absence (blue) and presence (red) of TCP— /3. The fact 
that all vectors point into the first quadrant is due to a minor technicality. B The variation of the 
squared radii is computed according to equation ^ for six different movies. The control cases are showrn 
in blue, the treatment cases in red. These figures are representative for three independent realisations of 

the experiment. 
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Figure 4. A Simulation of untreated cells growing according to the modified Fisher— Kolmogorov 
equation ^ in a square subfield of the culture dish with side 300 fim. Periodic boundary conditions are 
assumed to account for a balance of cells entering and exiting the subfield. Initially, 75 cells were seeded 
randomly across the domain. Snapshots are taken at t = 72h (initial datum, upper left) and t = 92, 112 /i 
and t = 132/i. The parameters are Dm = 5 fim'^h~^ , Dp = O.l fim^h~^ and a = 0.07/i~^. B Normalised 
growth curves (solid curves, simulation) and experimental measurements (discrete symbols) over a time 
course of 5 days. Shown are the untreated cells (control, +) and cells treated with TGF-/3 (o). 



